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interactions is investigated numerically. When the dot is occupied by two electrons we find that a 
level crossing between the two lowest many-body eigenstates may occur as a function of the spin- 
orbit coupling strength. This level crossing is accompanied by a non-vanishing magnetization of the 
ground-state. Using first order perturbation theory as well as exact numerical diagonalization of 
small clusters we show that the tendency of interactions to cause Stoner-like instability is enhanced 
by the SO coupling. The resulting g-factor can have a significant value, and thus may influence 
g-factor measurements. Finally we propose an experimental method by which the predicted phe- 
nomenon can be observed. 
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I. INTRODUCTION 

The effect of spin-orbit (SO) coupling on the energy 
spectrum of quantum dots (QDs) and metallic grains has 
attracted notable attention in the recent years^^. Much 
experimental and theoretical effort has concentrated on 
the magnetization of mesoscopic samples. For exam- 
ple, measurements of the g-factors of nano-particles using 
tunneling spectroscopy^! 6 - have led to several theoreti- 
cal studies 7 -^ which treated the electrons in the quan- 
tum dots as non-interacting particles. Other theoretical 
studies have considered interactions as well, while inves- 
tigating the interplay between interactions and disorder 
in quantum dots without SO coupling. It was shown 
that the combination of these effects can lead to non- 
trivial spin polarization ^ 10 ' 11 ! 12 ] 13 ! 14 ] 15 . In this paper we 
present an additional mechanism which can lead to spin 
polarization and non-trivial g-factor, where the role typ- 
ically played by disorder is taken by the SO coupling. 

Usually the g-factor is defined through the splitting 
of the Kramers doublet s 16 ! 17 in the presence of a weak 
magnetic field. Namely, the g-factor of the i-th single- 
particle level with spin a is given by 



where e ia (e\J ) is the corresponding energy level in the 
presence (absence) of a weak magnetic field H and fj, b is 
the Bohr magneton. The spin index a G {+, — } is used to 
denote the two time-reversed states according to the sign 
of the z component of their average magnetic moments. 
In the absence of a magnetic field each level is two-fold 
degenerate, and this degeneracy is lifted by the magnetic 
field, which increases the energy of one of the levels and 
decreases the energy of the other. Therefore, g^„ as de- 
fined by this formula can have either sign, depending on 
the direction of the energy change. The ground-state en- 
ergy always decreases when a magnetic field is applied, 
thus the g-factor of the ground-state obtained by Eq. (Q]) 



is positive. Usually, the value of g does not depend on 
the spin index, at least to zeroth order in iJ, so that one 
can denote the g-factor of the i-th level as ±<?j, with the 
convention that gt > 0. 

For free electrons the g-factor is constant, gi — 2 for 
each level i, and this value is more or less correct also 
for bulk measurements in various metals^-. However, in 
experiments performed on metallic nano-particles, values 
which are significantly less than the free value of the g- 
factor were obtained^. Moreover, large fluctuations in 
the measured values were seen. These findings attracted 
much theoretical attention, and resulted in studies which 
have obtained, within the framework of the random ma- 
trix theory (RMT), a description of the g-factor prob- 
ability distribution in the presence of SO coupling and 
disorder but in the absence of interactions 7 '^. In a re- 
cent work, the statistical properties of these distribution 
functions were related to several physical observable^. 
According to these results, the SO coupling influences 
the probability distribution of the g-factors of the dis- 
crete energy levels. The distribution function was shown 
to be universal, where the width is expressed in terms of 
various physical parameters. The presence of strong SO 
coupling and disorder results in sample to sample fluctu- 
ations of the g-factor. Moreover, the g-factor is expected 
to fluctuate also between different levels of a specific sam- 
ple, with a distribution function determined by RMT. 

Indeed, recent measurements of nano-particles have 
obtained g-factors which seem to be consistent with RMT 
predictions. For example, several experimental studies 
of metallic three dimensional nano-particles have shown 
the reduction of the measured g-factor as a function of 
the spin-orbit coupling strength. For aluminum nano- 
particles, in which the SO coupling is negligible, the 
measured g-factor values are approximately those of free 
electrons^ (i.e., g ps 2), while for gold nano-particles, in 
which the SO coupling is strong, the measured g-factors 
were in the range of 0.28— 0.45. 5 Furthermore, by extract- 
ing several g-factors from each sample, Petta and Ralph 
have succeeded to present an impressive confirmation of 
the theoretical RMT distribution function 6 . 
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Nevertheless, according to Eq. (Q} the g-factor mea- 
surement should compare the specific single-particle en- 
ergy level before and after the magnetic field is applied. 
However, practical experiments usually differ from that 
approach in two points. First, measurements are usu- 
ally related to the total energy of the system, and not to 
that of a specific level. Second, the measurement of the 
energy is sometimes indirect, as is the case in tunneling 
spectroscopy. 

These two points can be ignored if one neglects the 
interaction between particles. For a non-interacting sys- 
tem with an odd number of electrons, n e = 2p + 1, the 
change of the total ground-state energy due to magnetic 
field is equivalent to that of the highest filled level. n e — l 
electrons populate p Kramers pairs, where in each pair 
one level increases and the other decreases in the pres- 
ence of a magnetic field, so that their total contribution 
vanishes. The only contribution to the g-factor comes 
from the single electron occupying one level of pair p+l, 
so that if we define the g-factor of the ground-state with 
n e electrons by 



l(n e 



Ei°J(n e )-E£>(n e ) 



p, B H 



(2) 



where E gs (n e ) represents the total ground-state energy 
in the presence of a magnetic field H, then g(2p + 1) = 
g p +i- In addition, when the number of electrons n e is 
even, the total ground-state energy is not expected to 
change when a magnetic field is applied, since all the filled 
levels divide into pairs, in which the movement of one 
level is compensated by the other (to first order in H). 
Therefore, for an even number of electrons, a calculation 
of the g-factor using Eq. (f2]) gives g(2p) = 0. 

The second point, regarding the indirect energy mea- 
surement, requires an interpretation of the experimental 
results. For example, using tunneling spectroscopy one 
measures the gate voltage value for which a conductance 
peak of a QD occurs. At such an event the energies of the 
QD with n e — 1 and n e electrons and the gate voltage V g 
are related by the equation eV g = E gs {n e ) — E gs (n e — 1). 
When a magnetic field is applied, the position of the peak 
will change as a function of H. Therefore, by denoting 
the measured g-factor by g, one can analyze the peak 
motion in order to determine the g-factor, by calculating 
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2[eV g (0)-eV g (H)} 
/j, B H 



g{n e ) ~ g{n e -l). (3) 



Since either n e or n e — 1 is even, its corresponding g-factor 
vanishes, and thus g is equivalent to the other g-factor. 
Namely, g = g(n e ) or g — —g(n e — 1). Actually, since 
each peak is split in the presence of a magnetic field into 
two peaks having an opposite magnetic field dependence, 
extracting g from successive peaks results in the set of the 

single-particle g-factors, i.e. gi, -gi, g 2 , -52, 

As mentioned above, measurements done using tun- 
neling spectroscopy have indeed obtained g-factors which 



can be interpreted using RMT predictions. Nevertheless, 
as we have discussed, Eq. © is equivalent to Eq. (TT]) 
only for systems of non-interacting particles. Once inter- 
actions between electrons are important, it should be em- 
phasized that Eq. ([2]) is a definition of a many-particle 
g-factor, which depends on the total magnetization of 
the ground-state wave-function. For example, one can 
obtain spin contribution to the g-factor which is larger 
than 2, a phenomenon that cannot happen for a single- 
particle g- factor ^ i 10 i 1:L i 12 i 13 i 14 i 15 Indeed, by adding an in- 
teraction term to the RMT Hamiltonian, an increase of 
the g-factor fluctuations was reported) 19 ' 20 It was shown 
that the interactions result in a possibility of getting non- 
trivial spin values in the ground-state, and accordingly in 
an optional enhancement of the g-factor to values greater 
than 2. 

Although the theoretical studies of Refs. [l!^ and HO 
were performed for an odd-electron occupation, their re- 
sults suggest the possibility of a non-trivial spin polariza- 
tion for the even-electron case as well. If, for any reason 
the g-factor of an even-electron ground-state indeed dif- 
fers from zero, then the quantity measured in tunneling 
spectroscopy may not equal the single-level g-factor nor 
the many-particle g-factor. In such a case it should be 
related to the difference between two many-particle g- 
factors, as shown in Eq. 

In principle, the above description of the g-factor holds 
for both metallic and semiconducting dots. However, in 
semiconducting dots the strength of the SO coupling can 
be tuned by use of a gate voltag o 21 i 22 . Furthermore, sev- 
eral significant implications of SO effects in semiconduc- 
tors, such as spin-polarized field effect transistor— and 
spin Hall effec t 24 ! 25 , have recently attracted notable at- 
tention. 

With this in mind, we investigate in the current pa- 
per the ground-state magnetization properties of semi- 
conducting QDs with SO coupling where interactions 
between the electrons are considered. We show that 
the interplay between the SO coupling and the electron- 
electron interactions may result in a level crossing (LC) 
between the two lowest many-body levels. When these 
states are close in energy, the magnetic field splits them 
into two polarized states with a finite magnetization. As 
a result, there is a possibility to have a significant g-factor 
in the two-particle ground-state. Finally, we propose an 
experimental method which can be used in order to ob- 
serve the predicted phenomenon. 

We note that we have neglected so far in the intro- 
duction the orbital effect and its influence on the mag- 
netization of nano-particles. For three-dimensional (3D) 
nano-particles this is reasonable^. On the other hand, 
for two-dimensional (2D) systems the orbital effect is ex- 
pected to play an important role. For example, due to 
the orbital effect, the single-particle g-factor can exceed 
the value of 2. In addition, the g-factors of two levels be- 
longing to the same Kramers pair, gi t+ and gi.-, might 
be different. As a result, the g-factor of the doubly occu- 
pied ground state, which can be simply written, when the 
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electrons are non-interacting, as g(2) = gi. + + may 
not vanish. However, this contribution to g(2) which is 
the quantity of interest in this work, is linear in the mag- 
netic field H, and thus can be neglected for the weak 
fields used in such measurements. 

The rest of the paper is organized as follows. In the 
next section we describe the model Hamiltonian we use in 
order to incorporate, beside the magnetic field, both SO 
coupling and interactions between electrons. In section 
IIIII we present results for a system with non-interacting 
particles, which are shown to reproduce some known 
ground-state properties. In addition, wc find that there 
are specific values of the SO coupling strength, in which 
the Kramers doublet remains degenerate even when a 
magnetic field is applied. The effects of interplay between 
SO and electron-electron interactions are considered in 
section IIVI Our results point out that a finite magneti- 
zation can be obtained for systems with an even-particle 
occupancy. In section [V] we discuss the experimental rel- 
evance of this finding, i.e., the possibility that it might 
affect practical g-factor measurements. 

II. MODEL 

In order to model the semiconducting QD we use a 
tight-binding description of a finite 2D lattice with A 
columns and B rows (the number of sites is denoted by 
N = AB), with open boundary conditions, which is oc- 
cupied by n e electrons. As a result of a coupling between 
the spin degree of freedom and the orbital motion a finite 
probability for spin-flips during hopping processes exists. 
Separating the interactions from the free part, one can 
write the Hamiltonian as -ffQD — H Q + H- lnt , where the 
free part in the absence of disorder can be divided to a 
hopping term and a Zeeman term, i.e. Hq = H nop + Hb- 
The hopping part of the Hamiltonian is 

-ffhop = - ( V x a) mniJ a m ^ + i^, (4) 

+ Vy <Hn,n,<T®rn+l,n,<T' + H.C.), 

where n a (a m ,n,a) 1S a creation (annihilation) opera- 
tor of an electron with spin a in the lattice site placed in 
row to and column n. The matrices V x and V y are given 
by the Ando model 2 ^, which is the discrete version for 
the Rashba spin-orbit coupling 2 ^, as 

where V± (V2) is the hopping matrix element, for events 
which conserve (flip) the spin. The overall hopping am- 
plitude, t = \/Vx + V 2 2 , is taken as the energy unit of 
the problem. In other words, all energies are expressed 
in terms of t. 

The strength of the SO coupling can be expressed by 
the ratio between the absolute value of the spin-flip am- 
plitude and that of the total hopping element. Using a 



dimensionless parameter A = , 2 = Vilt. we exam- 

ine the entire range of A, between very weak (A — ► 0) 
and very strong (A < 1) spin-orbit coupling. Realistic 
values for A are between and 0.5i 26 ' 28 As mentioned 
above, these values can be controlled by tuning the gate 
voltage 2 ^. 

We now add a perpendicular magnetic field to our 2D 
sample, and we choose a gauge in which the vector poten- 
tial is A = —Hyx. The Zeeman term in the Hamiltonian 
is thus diagonal in spin space, and can be written as 

Hb = HbH ^ aa] n n (J a m ^ n ^, (6) 

m,n,er 

where a = ±1. 

With the gauge chosen, one has to modify the hop- 
ping elements in the x direction, according to the Peierls 
substitution^, and write V x — > V x e~ l0m . In this expres- 
sion m is the row number and 9 is a phase, that can be 
written as — 21T ^ S , where s is the lattice constant and 
0o = he/ e is the magnetic flux quantum. Thus, 9 is a di- 
mensionless parameter, that measures the magnetic flux 
throughout a lattice unit cell, in units of the quantum 
flux </>o. 

The Zeeman energy can be related to the hopping 
phase 9 and to the hopping amplitude t by the follow- 
ing considerations. One can express the absolute value 
of the Zeeman energy as \ibH — fJ.B<Po 2 ^ s ' 2 1 Substituting 

the physical constants Hb<P0 = 7^-1 where m is the elec- 

1 2 

tron mass, and using the relation t = tt- 2 — 1 , where m a 
is the effective mass at the bottom of the band, one gets 
I-IbH = = 9t 7 -^. The factor ^ depends on the 

' ^ 2mos^ mo mo 

specific type of the QD, and in general, m e ff « mo for 
metallic grains while m c g < mo for semiconducting ones. 
In the current study we use 2D geometry which is suitable 
for modeling typical semiconducting QDs. Moreover, the 
Ando model which incorporates the spin-orbit coupling 
was originally proposed for surfaces of III-V compound 
semiconductors 26 . Thus we set for the rest of this paper 
ss 4k < as in the case of GaAs. However, we have 

mo 15 ' 1 

checked that tuning this value does not lead to a qualita- 
tive change of the main results. Finally, since all energies 
are measured in units of t, the strength of the Zeeman 
term usH/t determines exactly the hopping phase. 
At last, the interactions term in the Hamiltonian is 

Hi n t = a m,n,T" m .» ! T"Ln,| am >"^' ( 7 ) 

m,n 

which represents a Hubbard interaction with strength U. 

The Hamiltonian Hqd is exactly diagonalized using 
the Lanczos procedure, for lattices of up to 15 x 14 sites, 
occupied by 1 or 2 electrons, and its lowest eigenstates 
are numerically found. In order to calculate the spin po- 
larization of the QD we apply a weak magnetic field along 
the z axis and calculate the expectation value of S z for 
the lowest levels. For the g-factor calculations, we com- 
pare the ground-state energies with and without the mag- 
netic field for each sample, and use Eq. @. The strength 
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of the magnetic field we apply is fisH/t ~ 10~ 4 — 10~ 3 , 
and for an experimental system in which the mean level 
spacing is 0.1 — 1 meV, it is equivalent to a magnetic field 
of 10 — 1000 G, in correspondence with realistic measure- 
ments. 



III. NON-INTERACTING ELECTRONS 

We start with non-interacting particles, by taking 
U = 0. Without the magnetic field, all single-particle 
states (and in particular the ground-state) are doubly- 
degenerate (the Kramers degeneracy) 16 ' 17 . When a mag- 
netic field is applied, it splits this degeneracy, and one 
gets to zeroth order in the magnetic field, (Si 1 ^) = 



— ('S'l 2 '), where {S z m> ) denotes the expectation value of 
the operator S z in the m-th eigenfunction (m = 1 being 
the ground-state). For A — > 0, \{S Z )\ -> \. When the SO 
coupling increases, a general decrease of \{S Z } \ can be ex- 
pected, and this trend can be seen in the upper panel of 
Fig.ffl 
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FIG. 1: (Color online) {Sz) of the lowest two single-particle 
levels (upper panel) and AE, the energy difference between 
them (lower panel, notice the semi- logarithmic scale), calcu- 
lated for a system of 8 x 7 sites in the presence of a magnetic 
field, as a function of the SO coupling strength. The value 
of {S z ) switches abruptly between the two levels (one level 
is shown by symbols and the other by a line) near A = 0.12, 
where the energy difference vanishes, implying a level crossing 
between the lowest two levels. The energy of the third level, 
however, remains much higher (lower panel, dashed line). 

However, one can see that the value of (S z ) switches 
abruptly between these two levels near A = 0.12. This is a 
sign of a level crossing (LC) , which can be seen by looking 
at the energy difference between these levels (lower panel 
of Fig. [1]). The switching of {S z ) occurs exactly when the 
energy difference vanishes. We note that such crossings 
occur also for large values of A. 

It is important to notice that the LC presented here 
occurs between states which belong to the same Kramers 



pair (i.e., states which are the time reversal of each 
other). The energy difference between states of differ- 
ent pairs is much larger, and although it is reduced by 
the SO coupling, yet it is usually much larger than the 
contribution of the weak magnetic field we apply to the 
energy (see the dashed line in the lower panel of Fig. [1} . 
As a result, crossings between states which originate from 
different pairs are much less probable. 

In the next section we study the case of doubly- 
occupied systems. For non-interacting electrons, based 
on the results of the current section, it is clear that a 
LC between the lowest two doubly-occupied states is im- 
probable, since their Slater determinants contain single- 
particle states from different Kramers pairs. Neverthe- 
less, as will be shown in the next section, the electron- 
electron interaction can change this picture qualitatively. 



IV. INTERPLAY BETWEEN INTERACTIONS 
AND SPIN-ORBIT COUPLING 

We now turn to study the effect of interactions on the 
behavior of the g-factor in the presence of SO. Calculat- 
ing the ground-state energies of the two lowest doubly- 
occupied many-body states, one finds that there is a LC 
between these states, at a certain value of the SO cou- 
pling (denoted in the following by A c ), as can be seen 
in the upper panel of Fig. O In the vicinity of A c , the 
expectation value of S* 2 switches smoothly between these 
states, as is shown in the lower panel of Fig. O 
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FIG. 2: (Color online) Typical results of the level crossing 
of the two lowest doubly-occupied states. The results shown 
were obtained for a system of 8 x 7 sites, with U = 3t. Up- 
per panel: the energy difference Ei — E\ is shown as a func- 
tion of the spin-orbit coupling strength A (notice the semi- 
logarithmic scale). The dip shows the crossing point. Lower 
panel: the switch of {S 2 ) between these two states, which is 
centered at the same place. 

As noted in the previous section, such a LC does 
not exist for non-interacting two-electron states since 
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it involves levels belonging to different single-particle 
Kramers pairs. Moreover, in the cases when there is a 
LC in the non-interacting system, i.e., between single- 
particle levels belonging to the same Kramers pair, the 
g-factor vanishes at the LC point (to zeroth order of H; 
it has however a linear magnetic field dependence from 
orbital effects). On the other hand, in the case of in- 
teracting electrons we find that both states involved in 
the LC have a significant magnetization of zeroth order 
in H. The magnetization properties, i.e., (S z ) and the 
g-factor, do not present a smooth switching as for (S 2 ) 
in the vicinity of the LC. Instead, when the energies of 
the two states become close enough to each other so that 
the energy associated with the magnetic field becomes 
important, both states develop a spin polarization as can 
be seen in Fig. [31 This leads to an enhancement of (S z ) 
in the crossing region, and to significant values of the 
g-factor. 
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FIG. 3: (Color online) Typical results of the spin polarization 
(S z ) (upper panel) and the g-factor, calculated using Eq. @ 
(lower panel), of the two lowest doubly-occupied states, in the 
regime of the level crossing between them. The results shown 
were obtained for a system of 8 x 7 sites, with U = 3t and 
H B H = 10- 4 t. 

The significant values obtained for the g-factor are cru- 
cially related to the degeneracy point (the LC) . Far from 
this point, when the two lowest many-particle states are 
not degenerate, each of these states tpi and ip2 is the 
time reversal of itself, i.e., T(ipx) = ip\ and T{i\)%) = V->2, 
where T is the time reversal operator. This immediately 
implies (S z ) — for both states, and the correspond- 
ing g-factors vanish as well. This picture changes in the 
vicinity of the degeneracy point, where the magnetic field 
breaks the degeneracy by polarizing both states. This of 
course results in a finite value of the g-factor. Such non- 
vanishing g-factor values can thus be seen as long as the 
energy associated with the magnetic field is larger than 
the level separation. Accordingly, as the magnetic field 
is enhanced, the peaks in (S z ) and g get wider. 

The dependence of the energy on the magnetic field 



is shown in Fig. [4l with a comparison between the LC 
regime to an arbitrary point. In the latter, a quadratic 
dependence of the ground-state energy on the magnetic 
field is clearly seen. On the other hand, near the LC point 
each of these states has a significant magnetization, and 
the dependence of the energy on the magnetic field is 
linear, with a finite value of the g-factor. 

These phenomena can be given a simple interpretation. 
Kinetic energy considerations make it advantageous to 
put the two electrons in the same orbital level, and cre- 
ate an unpolarized ground state. Repulsive interactions, 
however, cause a polarized ground state to be preferred, 
since the Pauli principle then tends to separate the elec- 
trons. Usually, the kinetic energy wins. However, SO 
coupling tends to reduce the single-particle level spacing 
so at some point the interactions win, and a Stoner-like 
instability emerge d 30 ' 31 . In order to support this intu- 
itive picture, we have calculated the energy difference 
between the two lowest many-body levels using first or- 
der perturbation theory in the interaction strength. As 
can be seen in Fig. [5J the interaction tends to decrease 
the energy difference between the lowest two many-body 
levels. When the electrons are non-interacting, the levels 
approach each other with increasing SO coupling, yet the 
minimal distance between them is much larger than the 
magnetic energy. The presence of interactions enhances 
this tendency, towards the situation in which a LC is 
possible. 
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FIG. 4: (Color online) The dependence of the energy on the 
magnetic field is compared between the regime of a level cross- 
ing (lower panel), to another arbitrary point (upper panel). 
The results shown by symbols were obtained for a system of 
8x7 sites, with U — 3t, and the solid lines represent quadratic 
(upper panel) and linear (lower panel) fits. 

From these results one can conclude that whereas the 
g-factor of a doubly-occupied system can be neglected for 
most values of A, it nevertheless has a significant value 
near A c . As is shown in Fig. when the system size 
increases the instability and the g-factor peak occurs for 
smaller values of the SO coupling. In the next section 
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we argue that such g-factor values might be significant 
even for realistic sample sizes and physical parameters, 
and thus they should not be neglected when analyzing 
experimental data. 
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FIG. 5: (Color online) The dependence of the energy dif- 
ference between the two lowest many-body levels on the SO 
strength for systems of 8 x 7 and 15 x 14 sites. The results 
shown are of the non- interacting case (dashed curve), and of 
interacting electrons with U = t calculated exactly (solid line) 
or by using first order perturbation theory (dashed-dotted 
line), notice the semi-logarithmic scale. 



V. EXPERIMENTAL RELEVANCE AND 
DISCUSSION 

In order to check whether the g-factor peak presented 
in the previous section occurs for realistic systems, one 
must study how the system size modifies this behavior. 
When the system is enlarged, one must be careful to leave 
the other physical parameters unchanged. The strength 
of the interactions is usually described by the parame- 
ter r s , which is defined through the ratio between the 
potential and the kinetic energies. The kinetic energy 
per electron in 2D samples goes like the electron den- 
sity n — n e /N. For Coulomb interactions one has E p ~ 
n jjUc/rdxdy per electron, where Uc is the Coulomb 
interaction strength between sites separated by one lat- 
tice constant. Since r ~ n -1 / 2 one gets E p ~ Uc\/n 
and thus r s ~ Ucl\pn. However, for Hubbard interac- 
tions E p ~ n J J U S(x — xo)S(y — yo) dx dy = nil, so that 
r s ~ U. Therefore, in order to keep r s constant, the 
value of U should stay unchanged when the system size 
increases. 

In Fig. [6] we show the dependence of the g-factor on A, 
for various system sizes, ranging from 8 x 7 to 15 x 14. 
As can be seen, the enhancement of the g-factor occurs 
for all of the curves, with some quantitative changes in 
the position and the height of the peak. Although a 
substantial enlargement of the system is not numerically 



possible because of the limitations of the exact diago- 
nalization technique, yet the trend is clearly seen. The 
value of A c which is found to decrease with increasing 
system size (see the inset of Fig. [5]) suggests that for 
a sufficiently large system size the crossing occurs for a 
moderate value of the SO coupling, which may be exper- 
imentally relevant. In addition, the modest increase of 
the peak height suggests that a significant peak may be 
observed for realistic system sizes. 
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FIG. 6: (Color online) The g-factor of the doubly-occupied 
ground-state as a function of the SO coupling strength A, for 
lattices sizes of 8 x 7, 10 x 9, 12 x 11 and 15 x 14 (size increases 
from right to left), with U = 3i. Inset: the dependence of A c 
on the number of sites. 

Finally, we would like to discuss the implications of 
the g-factor peak on g-factor measurements. Once a fi- 
nite magnetization of the doubly occupied ground state 
is possible, it can affect experiments done by, e.g., tun- 
neling spectroscopy. Such a measurement presents the 
result for the difference between two g-factors, as given 
by Eq. © . If the even-electron state has a non- vanishing 
g-factor, like in the vicinity of the LCs we have presented, 
the measured quantity g may not be equal to the g-factor 
of the state with an odd number of electrons, to which it 
is usually attributed. 

In such cases, a signature of the LC may be seen exper- 
imentally. In the regular case (as opposed to the LC sce- 
nario), the two levels which belong to the same Kramers 
doublet have the same g-factor up to a sign, and the 
dependence of the two energies on the magnetic field is 
symmetric. However, in the region of a LC, the two levels 
receive contributions from different even-particle states. 
Explicitly, with a magnetic field, the p-th Kramers pair is 
split to levels with different g-factors, g(2p— 1) — g(2p— 2) 
and g(2p) - g(2p - 1). Thus, if g(2p) or g(2p - 2) (or 
both) are not negligible, the magnetic field dependence 
will not be symmetric. Furthermore, since the strength 
of the spin-orbit coupling can be tuned by using a gate- 
voltag o 21 ' 22 , different shapes of the magnetic-field depen- 
dence may be obtained for a specific sample with differ- 
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ent values of the gate-voltage. An example is presented 
in Fig. [7] As one can see, the clearest non-symmetric 
behavior is obtained for A«A C (right panel), but such a 
dependence can be seen for a region in its vicinity as well 
(middle panel). Far enough from this region (left panel) 
the symmetric dependence reappears. 



h=0.51 M).54 h=0.57 




[E <0, -E (H, ]/t 

FIG. 7: (Color online) The magnetic field dependence of the 
first Coulomb peaks (n e = 1,2) for a lattice of 11 x 10 with 
U = 3i (for which A c « 0.59). 

We have also tried to verify that the reported phe- 
nomenon occurs for states with even electron numbers 
larger than 2 as well. However, the treatment of such 
cases is more difficult since the size of the Hilbert space, 
( 2 « B ) ' quickly passes the computational limit when n e 
increases. Therefore, the numerical simulation is limited 
to much smaller lattices, and although they show features 



similar to LC and enhanced g-factors reported for the two 
electron ground state, the question whether a LC occurs 
for larger lattices as well needs further investigation. As 
a possible method for that calculation we suggest the 
particle-hole version of the density matrix renormaliza- 
tion group algorithm^, which is suitable for such finite 
Fermi systems^. In addition, since as mentioned above 
the interplay between interactions and disorder can also 
result in a non-trivial spin polarization, the combination 
of both disorder and SO coupling with interactions can 
enhance this finding. These two points deserve a separate 
investigation. 

To conclude, we have shown that the combination of 
interactions and spin-orbit scattering can lead to a mag- 
netization of states having an even number of electrons. 
This effect was explained using first order perturbation 
theory by the tendency of interactions to drive a Stoncr 
instability, which is enhanced by the SO coupling. By 
studying the behavior when the system size increases, it 
seems that such a result may be experimentally observed 
even for realistic sizes of QDs. Therefore it might be rel- 
evant for understanding g-factor measurements. Based 
on our explanation of the results, we believe that similar 
phenomena might be observed in metallic nano-particles 
as well. However, a numerical investigation for 3D sys- 
tems is quite difficult. 
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